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' We investigate nuclear matter at finite temperature and density, including the formation of light 

I clusters up to the a particle (1 < A < 4). The novel feature of this work is to include the formation 

of clusters as well as their dissolution due to medium effects in a systematic way using two many- 
body theories: a microscopic quantum statistical (QS) approach and a generalized relativistic mean 

Q field (RMF) model. Nucleons and clusters are modified by medium effects. While the nucleon quasi- 
particle properties are determined within the RMF model from the scalar and vector self energies, 
QQ ■ the cluster binding energies are reduced due to Pauli blocking shifts calculated in the QS approach. 

Both approaches reproduce the limiting cases of nuclear statistical equilibrium (NSE) at low den- 
sities and cluster-free nuclear matter at high densities. The treatment of the cluster dissociation is 
based on the Mott effect due to Pauli blocking, implemented in slightly different ways in the QS 
and the generalized RMF approaches. This leads to somewhat different results in the intermediate 
density range of about 10""^ to 10~^ fm~^, which gives an estimate of the present accuracy of the 
theoretical predictions. We compare the numerical results of these models for cluster abundances 
i—{ • and thermodynamics in the region of medium excitation energies with temperatures T < 20 MeV 

(~| ' and baryon number densities from zero to a few times saturation density. The effect of cluster 

formation on the liquid-gas phase transition and on the density dependence of the symmetry energy 
is studied. It is demonstrated that the parabolic approximation for the asymmetry dependence of 
(N ■ the nuclear equation of state breaks down at low temperatures and at subsaturation densities due to 

1^ , cluster formation. Comparison is made with other theoretical approaches, in particular those, which 

■ are commonly used in astrophysical calculations. The results are relevant for heavy-ion collisions 

' and astrophysical applications. 

^ , PACS numbers: 21.65.Mn, 26.50.-l-x, 21.30.Fe, 25.75.-q, 97.60.Bw 
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I. INTRODUCTION 

The composition and the equation of state (EoS) of nuclear matter, possible phase transitions or condensates 
are widely discussed in several areas of nuclear physics. They attain an increasing importance in astrophysics and 
cosmology. Many of the systems under consideration have large charge asymmetries and cover a broad range of 
densities. Therefore, a deeper understanding of the composition and the thermodynamical properties of nuclear 
matter, and, in particular, the density dependence of the symmetry energy, from very low to super-saturation densities 
is of great interest. Let us mention three examples for systems where the knowledge of the symmetry energy in a 
broad range of densities is of crucial importance: (1) the surface structure of exotic nuclei with large neutron excess or 
new exotic collective modes, (2) the structure and composition of neutron stars from the ultra-dense core to the crust 
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at subsaturation densities with varying asymmetry due to /3-equilibrium, and (3) core collapse supernovae, where at 
high densities the symmetry energy determines the energy of the shock, and at low densities it affects the nuclear 
composition, neutrino interactions and aspects of nucleosynthesis. 

This wide spectrum of applications illustrates the importance of gathering reliable information on the EoS of nuclear 
matter, and, specificially, the symmetry energy, since predictions of its density dependence differ strongly for different 
theoretical approaches. Therefore, considerable efforts have been taken to obtain constraints from observational data 
on all aspects of the EoS. A promising source of information arises from the study of heavy-ion collisions, where 
transient states of very different densities can be investigated and the asymmetry can be varied to a certain extent 
by choosing the collision system. Central collisions at high energies yield large compressions and provide information 
from observations of nucleon flow [l| and particle production as shown for kaons, e.g., in 0, Recent reviews of 
studies of the symmetry energy in heavy- ion collisions were given in [1, [^, 0| . 

Recently, particular interest has been devoted to the properties of nuclear and neutron matter at very low densities, 
down to the limit of zero density. Below saturation density, correlations are expected to become important and 
nuclear matter can become inhomogeneous. An aspect of this is seen in fragmentation as a signature of the liquid-gas 
phase transition. At even smaller densities, down to one hundredth or one thousandth of saturation density and at 
moderate temperatures, few-body correlations remain important. This results from the fact, that at low densities the 
system can minimize its energy by forming light clusters such as deuterons, or particularly strongly bound a particles. 
With increasing density such clusters will dissolve due to the Pauli principle. Thus there exists an interesting detailed 
evolution of the correlations and the composition in nuclear matter depending on the density and temperature, which 
presents a challenge to a theoretical description. 

In the laboratory, low-density matter occurs in the outer regions of heavy nuclei, in halo nuclei, in expanding hot 
matter from heavy-ion reactions, in the envelopes of core- collapse supernovae, and also in recently discussed low- 
density isomeres such as the Hoyle state of ^^C Q- Experimental information on the composition and behavior of 
very low density nuclear matter was recently obtained by Kowalski et al. Q from the observation of the light particles 
emitted in low-energy heavy-ion collisions. It was found that the fraction of light particles is substantial at very low 
densities, implying a finite symmetry energy in the limit of zero density. 

Nonrelativistic Skyrme Hartree-Fock and relativistic mean-field calculations revealed a tight correlation between the 
density derivative of the neutron matter EoS near 2/3 of the saturation density and the neutron skin thickness Ai?„p 
of heavy nuclei . This observation translates into a correlation of the density dependence of the symmetry energy 

at saturation density with Ai?„p [1, [HI l12, 13, 14, i^lIBl- The scheduled PREX experiment at JLab is expected 
to provide a precise value of the neutron skin thickness of ^°®Pb from the observation of parity- violations in electron 
scattering [l7|, hence providing an independent constraint on the density dependence of the symmetry energy. 

At densities beyond nuclear saturation constraints on the EoS are expected in increasing quality from neutron star 
observables as, e.g., masses, mass-radius relations, gravitational binding energy and the cooling behavior of neutron 
stars. A recent review by Klahn et al. [l8| discusses these constraints in confrontation with those from heavy-ion 
collisions. It demonstates that at present it is far from trivial to obtain an EoS which is consistent with all these 
available observational data. The study of exotic structures like pasta phases in the inner crust of neutron stars gives 
information on very neutron rich matter around normal densities jl£, 2Qi, 2il [2^ . [23j . 

The actual composition of very low density matter is also relevant for the investigation of various stages in supernova 
explosions as pointed out in a number of recent publications. It is known to affect the effectiveness of the neutrino 
reheating of the shock wave [2^]. Also, the stellar core collapse is mostly determined during the dynamical plunge 
phase at densities between 10^^ and 10^^ g/cm^, where temperatures between 10^" and 10^^ K are reached [25l |26|. 
The sensitivity of the collapse dynamics on the properties of matter in this density regime could strongly influence the 
structure and composition of the proto-neutron star [27t and possible gravitational wave and neutrino signals emitted 
in various stages during and after core bounce. 

Only very few models for the EoS are applicable in actual supernova simulations. The reason is found in the 
required wide range of temperatures, densities and asymmetries, which are either not available, e.g. in tabular form, 
or not covered by the model. The most frequently used EoS are those of Lattimer and Swesty [2g| and of Shen, Toki 
et al. p9| . The former one is based on an extended liquid drop model for the nuclei embedded in a nucleon and a- 
particle gas, while the latter one has been developed in the framework of the Relativistic Mean Field (RMF) approach 
with non-linear meson self-couplings 30]. Shen et al. use the Thomas- Fermi approximation to describe heavy nuclei 
embedded in a gas. They consider a particles as a separate species, using an excluded volume prescription to model 
the dissolution of a particles at high densities which accounts for medium effects only in a very global way. The model 
neglects other light clusters such as deuterons, tritons, and helions ('^He). 

The problem of cluster formation in low-density nuclear matter has been addressed recently based on a virial 
expansion [sij to obtain the EoS for nuclear matter at all asymmetries including nucleons and a particles ('^He) j3^ 
and in later work also tritons and helions [ssj (see also 34 1 for a closely related formulation in an S-matrix approach 
and the quasiparticle gas model [1^). The virial coefficients in the Beth-Uhlenbeck approach [s^ are given by the 
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cluster bound state energies and scattering phase shifts. In Ref. [S^ these were taken directly from experiment, thus 
providing an exact limit for the EoS at very low densities, where the scattering is not yet influenced by medium 
effects. These results are believed to be reliable for densities up to rigat/lOOO and not too small temperatures. They 
provide a benchmark for other calculations. 

The present paper emphasizes that correlations, in-medium modifications of cluster properties and mean-field 
effects have to be considered simultaneously in the description of low-density nuclear matter, since all these affect the 
thermodynamical properties of the EoS. The occurence of clusters also changes the symmetry energy, since the cluster 
correlations depend on the asymmetry of the system. Here we will restrict ourselves to matter in thermodynamic 
equilibrium at temperatures T < 20 MeV and baryon number densities n < 0.2 fm~'^, where the quark substructure 
and excitations of internal degrees of freedom of nucleons (protons and neutrons) are not important and the nucleon- 
nuclcon interaction can be represented by an effective interaction potential. 

In this work, we explore two approaches to the problem: one is a quantum statistical (QS) formulation based on 
the thermodynamic Green function method [37| . This approach makes explicit use of an effective nucleon-nucleon 
interaction. It allows us to account for medium effects on the cluster properties. The second approach is a generalized 
relativistic mean-field (RMF) model, where the medium modified clusters are introduced as explicit degrees of freedom. 
The two methods have their strengths and deficiencies. The RMF method is reliable to determine the nucleon quasi- 
particle propagator in the medium, and this information is introduced into the QS model. The QS model, on the other 
hand, can determine the medium modifications of the clusters, such as the mass shift and the (momentum-dependent) 
Mott densities, where the clusters get dissolved. These are then introduced in parametrized form into the RMF 
method. 

The particle number density n,-(T, /Xp, /i„) of protons (t = p) or neutrons (r = n) depends on the temperature T 
and the chemical potentials fir- In the QS approach it is obtained from the single-particle spectral function, which 
can be expressed in terms of the self-energy. This is the main quantity to be evaluated. Considering the ladder 
approximation (btI . Ibs} . the formation of bound states is taken into account in a Bethe-Goldstone equation which 
in the low-density limit reduces to an effective Schrodinger equation. Effects of the medium can be included in a 
self-consistent way within the cluster mean-field approximation (for references see [s^, [s^ [13l)- The bound state 
energies are also modified due to Pauli blocking in the correlated medium. An extended discussion of the two-particle 
problem can be found in [4l| . This generalized Beth-Uhlenbeck formulation accounts for medium effects suppressing 
correlations at high densities. It allows to determine the second virial coefficient and also contains the Brueckner 
approach to matter near saturation density. The approach has been extended to three and four-particle bound states 
in Refs. [H, H^]. The medium dependent shift of the cluster binding energies has been investigated in [4l,|43|. We 
emphasize again that this quantum statistical approach avoids the introduction of semi-empirical concepts such as the 
excluded volume mechanism to mimic in-medium effects. However, since the quasi-particlc propagator is introduced 
from the outside (i.e. from the RMF model) the back effect of the clusters on the mean field is not included. 

The RMF model, on the other hand, takes this back reaction fully into account. On the other hand, only the 
bound state contributions of the clusters are included (from the QS model). Thus the continuum contributions are 
missing, which, as we shall see, leads to an overestimation of the two-particle correlations. Thus, by investigating 
both these models we also obtain an estimate of the remaining uncertainties of theoretical approaches. Also, in the 
present work, in both approaches, we have not yet included the contribution of heavier clusters, that should appear 
in the intermediate density range before the matter becomes homogenous again at densities near saturation. This 
problem will be treated in a later work. 

Extending the quasiparticle approach including the formation of light clusters enables us to describe the smooth 
transition from the low-density limit, where the nuclear statistical equilibrium (NSE) or the virial expansion are 
applicable, to the region of the saturation density where mean-field concepts have been successfully applied. None of 
the existing approaches to model the EoS gives satisfactory results in both regions simultaneously. More precisely, 
the EoS of Lattimer and Swesty as well as the EoS of Shen, Toki et al. [29i] fail to reproduce the NSE in the 
low-density limit, while the EoS in the virial expansion [3l|,[3l,|3l| ignores medium effects on light clusters and cannot 
describe the dissolution of clusters at high densities. 

The outline of this paper is as follows: In Sect.[lTlwe review the QS approach to the EoS and put particular emphasis 
on the calculation of the medium modifications of the clusters. As one result we obtain a density and temperature 
dependent modification of the binding energies of the clusters. In Sect. IIIII we introduce our generalized RMF model 
with light clusters as explicit degrees of freedom. We use a RMF model with density dependent meson-nucleon 
couplings [46[, which was used very successfully to describe nuclear structure in a wide region of the nuclear chart 
and has also been tested in heavy-ion collisions. The medium dependent masses of the clusters lead to a coupling of 
the nucleon and cluster dynamics. We also show how the thermodynamic quantities, such as free and internal energy, 
pressure and entropy, are obtained as functions of density and temperature. In Sect. lIVl we discuss the composition of 
nuclear matter and present the thermodynamical quantities for symmetric nuclear matter in both approaches. We also 
compare to the NSE model which gives the correct low density limit. Of particular interest is the phase transition from 
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(partially) clusterized to the homogeneous medium which is considered in Sect.|Vl In Sect. I VII we discuss specifically 
the symmetry energy of nuclear matter as a function of density and temperature, which is drastically changed at very 
low densities due to the cluster correlations. In Sect. IVIII we finally compare the results for the a cluster fraction 
in the QS and generalized RMF models with previous approaches and discuss the advantages of the present ones. 
We close with an outlook on further work, which should finally lead to an EoS that can be used in a wide range of 
problems, including nuclear structure, heavy-ion reactions, and supernovae simulations. Throughout the paper we 
use natural units where h = c — — ^■ 



II. QUANTUM STATISTICAL APPROACH TO THE EQUATION OF STATE 

A. Single-particle spectral function and quasiparticles 

Using the finite-temperature Green function formalism, a non-relativistic quantum statistical ap proach can be 
given to describe the equation of state of nuclear matter including the formation of bound states [3^, |4l| . It is most 
convenient to start with the nucleon number densities nr(T, jlp^ fin) as functions of temperature T and non-relativistic 
chemical potentials fir for protons (r = p) and neutrons (r = n) , respectively, 

nriT,ilp,iln) ^ ^^{a\ai)Sr,ri ^2 J J (1) 

where D, is the system volume, {1} — {fci, cri, ri} denotes the single-nucleon quantum numbers momentum, spin, and 
isospin. Summation over spin yields the factor 2 and 

/azM = (exp{/3[c^-ZAp-(A-^)Mn]}-(-l)^)"' (2) 

is the Fermi or Bose distribution function which depends on the inverse temperature (3 — 1/T. The non-relativistic 
chemical potential fl-r is related to the relativistic chemical potential fir by fir = fir + 'm-r with the nucleon mass mr- 
Instead of the isospin quantum number t we occasionally use the mass number A and the charge number Z . Both the 
distribution function and the spectral function 5*1(1, cj) depend on the temperature and the chemical potentials flp, 
fin, not given explicitly. We work with a grand canonical ensemble and have to invert Eq. ([T]) to write the chemical 
potentials as functions of the densities nn , ri n . For this EoS, expressions such as the Beth-Uhlenbcck formula and its 
generalizations have been derived [33, [sl, |41[ . 

We consider both the total number densities of protons and neutrons, n*°* and njij°*, and the temperature T as 
given parameters. Alternatively, the total baryon density n — -I- n*"* and the asymmetry of nuclear matter 
S = (njij°* — np°*)/n = 1 — 2Yp are used. Yp denotes the total proton fraction. In addition to the frozen equilibrium 
where n^°^ and n'°* are given, we assume homogeneity and isotropy in space. Thermodynamical stability is considered 
in Sects. |V] and IIVBI In a further development, allowing for weak interactions, /^-equilibrium may be considered, 
which is of interest for astrophysical applications. In that case the asymmetry S is uniquely determined for given n 
and T. 

The spectral function Si{l,uj) is related to the self-energy T,{l,z) according to 

o n = 2ImS](l,^-»0) 

' ' [c^-S(l)-ReE(l,a;)]2 + [ImE(l,w-iO)]2 ' ^' 

where the imaginary part has to be taken for a small negative imaginary part in the frequency lu. E{1) — k\l(2mi) 
is the kinetic energy of the free nucleon. The solution of the relation 

£;5»"(1) =£;(!)+ Re I][l,£;5»"(l)] (4) 

defines the single-nucleon quasiparticle energies £^^"(1) = E{1) + Ai?^^(l). Expanding for small Im 5](1, z), the spec- 
tral function yields a J-like contribution. The densities are calculated from Fermi distributions with the quasiparticle 
energies so that 

n^"(r,A^p,Mn) = ^5]/i,z[i?r(l)] (5) 

ki 

follows for the EoS in mean field approximation. This result does not contain the contribution of bound states and 
therefore fails to be correct in the low-temperature, low-density limit where the NSE describes the nuclear matter 
EoS. 
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As shown in Refs. [3^, |41|, the bound state contributions are obtained from the poles of Im E(l, z) which cannot 
be neglected in expanding the spectral function with respect to Im z). A cluster decomposition of the self-energy 
has been proposed, see [3al- The self-energy is expressed in terms of the A-particle Green functions which read in 
bilinear expansion 

GAil...A, 1' . . . A', za) = ■ • ■ ^ Ka-ET (K) "^*^"''^^' ■ • ■ ^'^ ■ 



The A-particle wave function ipAvKi^ . . . A) and the corresponding eigenvalues E'^^{K) result from solving the in- 
medium Schrodinger equation (see the following subsections). K denotes the center of mass momentum of the 
j4-nucleon system. Besides the bound states, the summation over the internal quantum states v includes also the 
scattering states. 

The evaluation of the equation of state in the low-density limit is straightforward. Considering only the bound-state 
contributions, we obtain the result 

n;°*(T,/ip,/i„) = ^ ^ ZfAAEXu{K;Tr^ip,nn)], 

<°*(T,/2p,/i„) = ^ E {A~ Z)fAAETAK;T,lipr^n)] (7) 

A,v,K 

for the EoS describing a mixture of components (cluster quasiparticles) obeying Fermi or Bose statistics. The total 
baryon density results as n{T, flp, fin) — ^^^^{T, fip, jln) + «p°'(^7 Mp; An)- To derive the extended Beth-Uhlenbeck 
formula, see [42| . we restrict the summation to j4 < 2, but extend the summation over the internal quantum numbers i/, 
not only to the excited states, but also the scattering states. Note that at low temperatures Bose-Einstein condensation 
may occur. 

The NSE is obtained in the low-density limit if the in- medium energies E'Xj^{K;T, (xp, jln) can be replaced by the 

binding energies of the isolated nuclei e''°1{K) = + K^/{2Am), with m = 939 MeV the average nucleon mass. 
For the cluster contributions, i.e. A > I, the summation over the internal quantum numbers is again restricted to the 
bound states only. We have 



bound 



A,u,K 
bound 



nfE(T,Ap,M„) = ^ E (8) 

A,u,K 

The summation over A includes also the contribution of free nuclcons, A ^ 1. 

In the nondegenerate and nonrelativistic case assuming a Maxwell-Boltzmann distribution, the summation over the 
momenta K can be performed analytically and the thermal wavelength A — ■\/27r/ [mT) of the nucleon enters. As 
shown below, the medium effects in nuclear matter are negligible below 10^* times the saturation density rigat for the 
temperatures considered here. 

Interesting quantities are the mass fractions 

^^■^ ^ E fAMEfJK; T, fip, A„)] (9) 

of the different clusters. From the EoS considered here, thermodynamical potentials can be obtained by integration, 
in particular the free energy per volume F/fl. In the special case of symmetric nuclear matter, Yp = 0.5, the free 
energy per volume is obtained from the averaged chemical potential fl = {flp + /i„)/2 as 

F(T, n, Y;)/n - / dn' ^(T, n', Y^) . (10) 

In the quantum statistical approach described above, we relate the EoS to properties of the correlation functions, 
in particular to the peaks occurring in the A-nucleon spectral function describing the single-nucleon quasiparticle 
{A = 1) as well as the nuclear quasiparticles {A > 2). Different approaches to these quasiparticle energies can be 
given by calculating self-energies which reproduce known properties of the nucleonic system. In the following Subsects. 
Ill Bl and III CI we discuss results obtained from a microscopic Hamiltonian approach to nuclear matter. In Sect. IIIll a 
relativistic mean field (RMF) approach is given, which is based on an effective nucleon-meson Lagrangian. 
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B. Medium modification of single nucleon properties 

The single-particle spectral function contains the single-nucleon quasiparticle contribution, = E!}^{k), given 

in Eq. where r denotes isospin of particle 1 and k is the momentum. In the effective mass approximation, the 
single-nucleon quasiparticle dispersion relation reads 

Erik) = A£;f (0) + ^ + 0{k^) , (11) 

where the quasiparticle energies are shifted at zero momentum k by AE^^{0), and m* denotes the effective mass of 
neutrons (r = n) or protons (t = p). Both quantities, AEf^{0) and m*, are functions of T, Up and n„, characterizing 
the surrounding matter. 

Expressions for the single-nucleon quasiparticle energy Er{k) can be given by the Skyrme parametrization [47| or 
by more sophisticated approaches such as relativistic mean-field approaches [3^, see Sec. IIIH and relativistic Dirac- 
Brueckner Hartree-Fock [4^ calculations. We will use the density-dependent relativistic mean field approach of fi^ 
that is designed not only to reproduce known properties of nuclei, but also agrees with microscopic calculations in 
the low density region. It is expected that this approach gives at present an optimal fit to the quasiparticle energies 
and is applicable in a large interval of densities and temperatures. 

Microscopic calculations are based on a model describing the interaction between the nucleons. To go beyond the 
mean- field approximation, strong interaction as well as bound state formation has to be taken into account. This can 
be done in the low-density region where in the non-relativistic case a T-matrix can be introduced. We start from a 
nonrelativistic Hamiltonian in fermion second quantization 

H = "^E{l)a\ai + ^ ^ V{l2,l'2')a\ala2>ay , (12) 

1 12,1'2' 

where the kinetic energy is E{1) — / (2toi), and the potential energy contains the matrix element V{\2, 1'2') of the 
nucleon-nucleon interaction. 

Since there is no fundamental expression for the nucleon-nucleon interaction, a phenomenological form is assumed 
to reproduce empirical data such as the nucleon scattering phase shifts. Different parametrizations are in use. For 
calculations one can use potentials such as PARIS and BONN or their separable representations [io']. To obtain the 
empirical parameter values of nuclear matter at saturation density, three-body forces have been introduced in the 
Hamiltonian (|12p. In particular, the Argonne AV18/UIX potential [H^l has been used to calculate light nuclei [5l| . 

Replacing the two-particle T-matrix in Born approximation with the interaction potential V ^ we obtain the Hartree- 
Fock approximation for the energy shift 

Ai?HF(l) = ^[^(12, 12) - 1/(12, 21)]/i,.j£;(2) + Ai?HF(2)] . (13) 

2 

In this approximation, all correlations in the medium are neglected. The self-energy does not depend on frequency, 
i.e. it is instantaneous in time, with vanishing imaginary part. 

A full Dirac-Brueckner Hartree-Fock (DBHF) calculation has been performed by Fuchs ^] and has been compared 
with RMF approaches. The relation between the T matrix approach and the Brueckner G matrix approach was 
discussed in detail in Ref. [53| . Extended work has been performed using sophisticated interaction potentials to 
evaluate the quasiparticle energies in the DBHF approximation, for recent reviews see Refs. [ 48; . 52, 54]. There was 
reasonable agreement between the RMF parametrisation of the quasiparticle energies and the DBHF results. 

We can assume [l8| that the density-dependent RMF parametrisation covers a large density region (which will be 
discussed in detail in Sect. IIII[) and that it can be used instead of the above Hartree-Fock shifts to determine the 
single-nucleon quasiparticle energies. They result as 

E^piO) = ^[m - 5]„,p(T, n, 5)Y + P + S"„^p(r, n, 5) , (14) 

where and ^ are the scalar and the time component of the vector self energy, respectively. In the nonrelativistic 
limit, the shifts of the quasiparticle energies are 

AEllik) = E°,p(T, n, 6) - E„,p(T, n, 6) . (15) 

The effective masses for neutrons and protons are given by 



"^n,p ^m - E„,p(T, n, 5) 



(16) 
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Approximations for the functions S° _p(T, n, S) and E„_p(T, n, S) are given in the Appendix. These functions reproduce 
the empirical values for the saturation density nsat ~ 0.15 fm~'^ and the binding energy per nucleon B/A « —16 MeV, 
see Subsection IIII Dl The effective mass is somewhat smaller than the empirical value m* w m(l — 0.17 n/nsat) for 
n < 0.2 fm^3. 



C. Medium modification of cluster properties 



Recent progress of the description of clusters in low density nuclear matter [27|, l55|, l56|, l57[ enables us to evaluate 
the properties of deuterons, tritons, helions and helium nuclei in a non-relativistic microscopic approach, taking the 
influence of the medium into account. 

In addition to the i5-like nucleon quasiparticle contribution, also the contribution of the bound and scattering states 
can be included in the single-nucleon spectral function by analyzing the imaginary part of I](l,z). Within a cluster 
decomposition, vl-nucleon T matrices appear in a many-particle approach. These T matrices describe the propagation 
of the ^-nucleon cluster in nuclear matter. In this way, bound states contribute to rir = nr{T, fin, flp), see |38l. l4l|. 
Restricting the cluster decomposition only to the contribution of two-particle correlations, we obtain the so-called 
T2G approximation. In this approximation, the Beth-Uhlenbeck formula is obtained for the EoS, as shown in [H, |4l[. 
In the low-density limit, the propagation of the ^-nucleon cluster is determined by the energy eigenvalues of the 
corresponding nucleus, and the simple EoS ([7]) results describing the nuclear statistica l eq uilibrium (NSE). 

For nuclei imbedded in nuclear matter, an effective wave equation can be derived j38l. . The ^-particle wave 
function iPauk{^ . . . A) and the corresponding eigenvalues E'^^{K) follow from solving the in-medium Schrodinger 
equation 

[£;^"(1) + . . . + - EXAKMaM^ ...A) 

+ E Y.\^-hi)-fmV{^L^'J') n 5uk'^AuK{l' ■■■A')^Q. (17) 

V...A' Kj k^i,] 

This equation contains the effects of the medium in the single-nucleon quasiparticle shifts as well as in the Pauli 
blocking terms. The A-particle wave function and energy depend on the total momentum K relative to the medium. 

The in-medium Fermi distribution function /(I) = (exp{/? [£^^^"(1) — jli]} -I- 1) ^ contains the non-relativistic ef- 
fective chemical potential fli which is determined by the total proton or neutron densities (i.e. including those bound 
in clusters) calculated in quasiparticle approximation, = f^^^ X^i /(1)'^ti,t fo'" particles inside the volume 
fi. It describes the occupation of the phase space neglecting any correlations in the medium. The solution of the 
in-medium Schrodinger equation (|17p can be obtained in the low-density region by perturbation theory. In particular, 
the quasiparticle energy of the A-nucleon cluster with Z protons in the ground state follows as 

ETJK) = EXziK) = E% + ^ + ^EfziK) + AE^^^fiK) + AE^:^\K) + ... (18) 

with various contributions. Besides the cluster binding energy in the vacuum i?^^^ and the kinetic term, the self- 
energy shift AE^ziK)^ the Pauh shift AE^'^^f'iK) and the the Coulomb shift AE1°^\K) enter. The latter can be 
evaluated for dense matter in the Wigner-Seitz approximation [i^, HI, H^. It is given by 







f2np\ 
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"sat 



(19) 



with vq — 1.2 fm. Since the values of Z are small, this contribution is small as well and disregarded here in the 
quasiparticle energy p^ . 

The self-energy contribution to the quasiparticle shift is determined by the contribution of the single-nucleon shift 

A£^l^z(0) ^{A- Z)AEf^{0) + ZA£;SE(0) + ^£;SE^cff.ma.s ^20) 

The contribution to the self-energy shift due to the change of the effective nucleon mass can be calculated from 
perturbation theory using the unperturbed wave function of the clusters, see [l^j, so that 



AE^ z ^ ^ ^"^^^ ■ ^ ' 
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Values of sa,z for {A,Z} = {i} = {d,t,h,a} are given in Tab. HI Inserting the medium-dependent quasiparticle 
energies in the distribution functions ^ the first two contributions to the quasiparticle shift in (PO)) can be included 
renormalizing the chemical potentials. 

The most important effect in the calculation of the abundances of light elements comes from the Pauli blocking 
terms in Eq. (jl7[) in connection with the interaction potential. This contribution is restricted only to the bound states 
so that it may lead to the dissolution of the nuclei if the density of nuclear matter increases. The corresponding 
shift AE^^^^^{K) can be evaluated in perturbation theory provided the interaction potential and the ground state 

wave function are known. After angular averaging where in the Fermi functions the mixed scalar product k ■ K 
between the total momentum K and the remaining Jacobian coordinates k is neglected, the Pauli blocking shift can 
be approximated as 

AE^f'iK) « AiJ^ff (0) exp (~^^^) • (22) 

Avoiding angular averaging, the full solution gives the result up to the order 

AE^X:f{K) « AE^A^fiO) exp (-^) (23) 

with the dispersion that can be calculated from 

g,{T,n,Yp) = — — . 24) 

The values for gi^i and gi^2 can be calculated from perturbation theory using the unperturbed cluster wave functions; 
the density corrections hi^i and hi^2 are fitted to variational solutions of the in- medium wave equation Eq. (|17p for 
given T, n^, n„ and K. Numerical values of the parameters in symmetric nuclear matter {Yp = 0.5) are given in Tab 
III 

The shift of the bindin g en ergy of light clusters at zero total momentum which is of first order in density [ssl . [56j has 
been calculated recently [53|. The light clusters of the deuteron (d — ^H), triton {t = ■^H), helion {h = ^He) and the 
a particle (*He) have been considered. The interaction potential and the nucleonic wave function of the few-nucleon 
system have been fitted to the binding energies and the rms radii of the corresponding nuclei. 

With the neutron number Nj = Aj — Z, , it can be written as 



AEl^%{S)-np,nn,T) = -A [Z^n^ + N,n^]5E^-'''\T,n) , (25) 



where the temperature dependence and higher density corrections are contained in the functions 5Ej^^^^ (T, n) . These 
functions have been obtained with different approximations for the wave function. In case of the deuteron, the Jastrow 
approach leads to a functional form 



- ^/na^^3 exp (af ^yi) erfc (oi^sy^) 



1 + [bi,i + 6i.2/T] 



^ (26) 



with Hi = 1 + a.;. 2/7". For the other clusters i = t,h, a, the Gaussian approach is used which gives the simple form 

The parameters a^^i, ai^2 and 0,1^3 are determined by low-density perturbation theory from the unperturbed cluster 
wave functions. The parameters bi^i and bi^2 are density corrections and are fitted to the numerical solution of the 
in-medium wave equation Eq. ()17|) for given T, Tip , Tin^ P — . Values are given in Tab III 

Now, the nucleon number densities (jT]) can be evaluated as in the non-interacting case, with the only difference that 
the number densities of the particles are calculated with the quasiparticle energies. In the light cluster-quasiparticle 
approximation, the total densities of neutrons 

<°*=^«+ ^^'^^ (28) 

2— d,t,/i,a 
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TABLE L Parameters for the cluster binding energy shifts. 



cluster i 


Si 




ai,2 


a,i,3 bi^i 


bi,2 




gi,2 


hi.i 


hi,2 




[MeV] 


[MeV^/^fm^] 


[MeV] 


[fm^] 


[MeV fm^] 


[fm-^] 


[MeV-^fm'2] 


[fm] 


[fm^] 


d 


11.147 


38386.4 


22.5204 


0.2223 1.048 


285.7 


0.85 


0.223 


132 


17.5 


t 


24.575 


69516.2 


7.49232 


4.414 


43.90 


3.20 


0.450 


37 




h 


20.075 


58442.5 


6.07718 


4.414 


43.90 


2.638 


0.434 


43 




a 


49.868 


164371 


10.6701 






8.236 


0.772 


50 





and of protons 

nl°'=np+ J2 (29) 



tot 

= 7/„v^ -I- 

i—d^tM a 



contain the densities of the free neutrons and protons n„ and rip, respectively, and the contributions from the nucleons 
bound in the clusters with densities n^. The state of the system in chemical equilibrium is completely determined 
by specifying the total nucleon density n = + 7i*°', the asymmetry S and the temperature T as long as no 
/3-equilibrium is considered. 

This result is an improvement of the NSE and allows for the smooth transition from the low-density limit up to 
the region of saturation density. The bound state contributions to the EoS are fading with increasing density because 
they move as resonances into the continuum of scattering states. This improved NSE, however, does not contain the 
contribution of scattering states explicitly. For the treatment of continuum states in the two-nucleon case, as well as 
the evaluation of the second virial coefficient, see [3^. l4ll]. 

The account of scattering states needs further consideration. Investigations on the two-particle level have been 
performed and extensively discussed [H, [H, [4l[ . We use the Levinson theorem to take the contribution of scattering 
states into account in the lowest-order approximation. Each bound state contribution to the density has to accom- 
panied with a continuum contribution that partly compensates the strength of the bound state correlations. As a 
consequence, the total proton and neutron densities are given by 

^ bound 

n;°*(T,/ip,/i„) ^ t:J2 z[fAMEX.(K;T,fipr^n)]-fAMET::.\K-,T,fip,fin)]\ , (30) 

A,iy,K 
^ bound 

nl°\T,fip,lin) = o E (^-^) [,/Az[i?X.(i^;r,/ip,An)]-/Az[i?i:;i*(if;r,Mp,M„)]] (3i) 

Am,K 

with explicit bound and scattering terms. -E™^* denotes the edge of the continuum states that is also determined by 
the single-nucleon self-energy shifts. These expressions guarantee a smooth behavior when the bound states merge 
with the continuum of scattering states. The summation over A includes also the contribution of free nucleons, A — 1, 
considered as quasiparticles with the energy dispersion given by the RMF approach. 

The summation over K and the subtraction of the continuum contribution is extended only over that region of 
momentum space where bound states exist. The disappearance of the bound states is caused by the Pauli blocking 
term; the self-energy contributions to the quasiparticle shifts act on bound as well as on scattering states. Above the 
so-called Mott density, where the bound states at i^T = disappear, the momentum summation has to be extended 
only over that region K > K^'^*-{T,n,6) where the bound state energy is lower than the continuum of scattering 
states. The contribution of scattering states is necessary to obtain the second virial coefficient according to the Beth- 
Uhlenbeck equation, see [3l|, |4l[ . This leads also to corrections in comparison with the NSE that accounts only for the 
bound state contributions, neglecting all effects of scattering states. These corrections become important at increasing 
temperatures for weakly bound clusters. Thus, the corrections which lead to the correct second virial coefficient are 
of importance for the deuteron system, when the temperature is comparable or large compared with the binding 
energy per nucleon. In the calculations for the quantum statistical (QS) model shown below, the contributions of 
these continuum correlations have been taken into account. 

Solving Eqs. ([50)1 and ([5T|) for given T, n^"* and 7i'°* we find the chemical potentials /Xp and After integration, 
see Eq. (jlOp . the free energy is obtained, and all the other thermodynamic functions are derived from this quantity 
without any contradictions. Results are given below. 

We do not consider the formation of heavy clusters here. This limits the parameter range n^°*, n*"', T in the phase 
diagram to that area where the abundances of heavier clusters are small. For a more general approach to the EoS 
which takes also the contribution of heavier cluster into account, see (45| . Future work will include the contribution 
of the heavier clusters. 
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Further approximations refer to the linear dependence on density of the shifts of binding energies, calculated in 
perturbation theory. A better treatment will improve these shifts, but it can be shown that the changes are small. 
The approximation of the uncorrelated medium can be improved considering the cluster mean-field approximation 
[ssl [sol . [57} . Furthermore, the formation of quantum condensates will give additional contributions to the EoS. 
However, in the region considered here the formation of quantum condensates does not appear. This is in contrast 
to a recent work employing a quasiparticle gas model [s^ where Bose-Einstein condensation of deuterons is observed 
because the Pauli shift of the deuteron binding energy at high densities is not considered. 

III. GENERALIZED RELATIVISTIC MEAN-FIELD MODEL WITH LIGHT CLUSTERS 

A main ingredient to construct the low-density EoS is the proper determination of the nucleonic quasiparticle 
energies that enter the single-nucleon distribution functions, Eqs. (|5l7p but also the cluster energies via the in- 
medium Schrodinger equation (jl7p . Recently, realistic values for the nucleon quasiparticle shifts were obtained from 
sophisticated calculations within Hamiltonian approaches, such as Dirac-Brueckner Hartree-Fock calculations fs^^'wl. 
RMF approaches proved to be very successful to interpret properties near saturation density, see, e.g., Refs. jl8i. |46|. 
We extract the single-nucleon quasiparticle shifts from the results of the RMF model with density-dependent couplings 
and use them in our QS approach. In the following we show how this RMF model can be extended to include light 
clusters, which are considered as quasiparticles modified by medium effects as obtained in the QS approach. A 
comparison of the generalized RMF model with the QS model will show distinct differences in the thermodynamical 
properties that are related to the employed approximations. 

In a conventional relativistic mean-field description [30l | of homogeneous and isotropic nuclear matter, nucleons 
interact by the exchange of mesons where usually isoscalar lo and <j and isovector p mesons are included. Neutrons 
and protons are described by Dirac spinors ■0^ (i = n,p). The mesons are represented by Lorentz vector fields ui^ 
and and Lorentz scalar fields a. The electromagnetic interaction is not considered in nuclear matter. A possible 
isovector, Lorentz scalar 5 meson is not included in the present model. The mesons couple minimally to the nucleons. 
In our approach, non-linear meson self-interactions are not introduced, but the couplings are assumed to be functionals 
of the nucleon operator- valued currents in order to simulate a medium dependence of the interaction. 

In the generalized RMF model with light clusters, the ground states of the deuteron (d = ^H), the triton (< = -^11), 
the helion {h — "^He) and the a particle CHe) are introduced as additional degrees of freedom with the corresponding 
spin field (pa, spin 1 field cff^ and spin 1/2 fields tpi {i = t,h). The clusters are treated as point- like particles and 
their internal structure is not taken into account. The influence of the medium on the cluster properties is described 
by density and temperature dependent shifts of the binding energies as introduced in the previous section. The Pauli 
shifts, cf. Eq. ([^5]) . are taken from the nonrelativistic calculation neglecting the dependence on the cm. momentum 
K of the cluster, while this is taken into account in the QS approach. On the other hand, the self-energy shift is 
treated self-consistently in the RMF description in contrast to the QS approach where it enters in parametrized form 
from an independent model, namely the RMF model decribed in this section. 

A. Lagrangian density and field equations 

In the present approach, the model Lagrangian has the form 

= J2 ^^{l^.^D^ ~M,)^, + ^{iD^ipa)*{iDa^,ipa)-^^lMlipa (32) 



with the field tensors 



G^i, = d^Uy - d^u^ Hf^^ = di_,p^ - d^Pf_, (33) 



of the Lorentz vector fields. Vectors in isospin space carry an arrow. Nucleons form an isospin doublet with r^tpn — V'n 
and T^^p — —ipp- Similarly, for the triton and helion one has r^ipt — 4't and Tsiph = —i^h, respectively. Deuterons 
and a-particles are treated as isospin singlets. 
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The covariant derivative 

= - T^AiUJ^" - Tp\Ni - Zi\T ■ (34) 

for a particle i contains the interaction with the Lorentz vector mesons with a strength that is determined by the 
density-dependent couphngs T^, Tp and the mass (A^), neutron (Ni) and proton {Z^) number of a particle i. The 
scalar a meson with coupling strength r^. appears in the effective mass 

Mi = mi- V^AiO - AB, (35) 

of a particle i with vacuum rest mass rrij. The vacuum rest mass of a cluster i = d,t,h,a is given by 

mi = ZiTUp + NiTUn - Bf , (36) 

which defines the vacuum binding energies B^ > 0. The medium dependent Pauli shift ABi appears only for clusters. 
The couplings = ^m{Q) {m = u),a, p) are functionals of the Lorentz scalar density 

Q = ^fim (37) 

that contains the free nucleon current 



J''=j^+jii (38) 

with = xpi'y^^xlji. The Pauli shifts of the binding energies were derived in the previous section as a function 
ASj (n*°* , n^°* , T) depending on the total proton and neutron densities n*"* and n^°* and the temperature T. In 
principle, the densities have to be replaced by the corresponding quantities expressed in terms of the field operators of 
the nucleons and clusters. In the case of the fermions, this poses no problem since the currents of the triton and helion 
have the same form as the currents of the nucleons. However, for the bosons, the definition of the conserved currents 
contains the meson fields and the Pauli shifts. A dependence of the Pauli shifts on these currents leads to additional 
rearrangement contributions that modify the self-energies and requires a redefinition of the conserved currents. In 
order to avoid these complications, the dependence of the Pauli shifts on the densities is replaced by a dependence 
on the vector meson fields. In this way the Pauli shifts are treated in an equivalent way as the usual single-particle 
shifts. In consequence, we replace the dependence on the total densities 

^ <=\ [Bo. + Qp] (39) 
by pseudo-densities and with the quantities 

= \>^^/uJ^^p (41) 

and 

ep = \\fW^>.- (42) 

The coefficients A^^ = m^,/rj^(0) and = m^/rp(0) are defined such that the correct low-density limit is obtained. 

The field equations for mesons, nucleons and clusters are derived from the Lagrangian density in the usual way. 
They are solved self-consistently in the mean-field approximation where the meson fields are treated as classical 
fields and sea-states of the fermions are not considered. The couplings Ym {rn, = lu, a, p) become simple functions 
of g = \/ {J^){Jp) where the brackets (•) indicate the summation over all occupied states of the system. The field 
equations simplify considerably due to the symmetries of homogeneous and isotropic nuclear matter at rest. The 
scalar meson field is directly given by 



with the source density 



(43) 

mi 
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that is a sum of the scalar densities nf = (tpifpi) of the fermions i = n,p,t,h and the scalar densities = {(p'^Mdifidn) 
of the deuteron and — {(paMa^a) of the a particle. The non- vanishing components of the vector meson fields are 

with two different source contributions. The regular source densities 

n„ = A,n, (47) 

i—n,p,d,t,h.a 

= E (7V,-Z,)n, (48) 

i—n,p,d^t,h^OL 

depend on the vector densities iii = (ipijo'ipi) of the fermions i = n,p, t, h, the vector density of the deuteron 

"-d = ^((j^doiySdM - iDf.ipdo)* Vd + Vd* {iDdofdti - iDdf^fdo)) (49) 

and of the a particle 

"a = ^{(iDaO'^a)* + [iDaQ^Pa)) ■ (50) 

The second contribution in eqs. (|45p and (|46|) with the derivatives of the binding energy shifts is proportional to the 
scalar densities of the clusters. 

The Dirac equation for the spin 1/2 particles {i — n,p,t, h) assumes the form 

[7*^ {id^ - S,^) - (m, - E,)] V, = (51) 

with scalar and vector self-energies Si and S^^, respectively. The field equations for the a-particle and the deuteron 
are the Klein-Gordon equation 

- {idf, - (z9^ - E„^) + (m„ - S„)'] = (52) 

and the Proca equation 

- - S,^) [{zd'^ - E^) - - E^) 'P'^d] + i^d " ^df Vd = ^, (53) 

respectively. The scalar self-energies are given by 

E, = r^A.fj + AS, (54) 

where the binding energy shift AB^ appears only for clusters. The non- vanishing component of the vector self-energy 
is the zero-component 

E„ = r^A,c^o + Tp (TV, - Zi) (po)3 + Sfo (55) 

with the 'rearrangement' contribution 

^io = C^o"^^ + (po)3 - r>7i^ (56) 

that appears only for nucleons. It contains contributions with derivatives = dTm/dg of the meson-nucleon 
couplings. 

Since the self-energies are momentum-independent in homogeneous and isotropic nuclear matter, the field equations 
of the nucleons and of the clusters are easily solved. The solutions are plane waves with shifted masses and energies as 
compared to the vacuum solution, i.e. nucleons and clusters can be considered as quasiparticles. At finite temperatures 
T the relevant vector and scalar densities are easily calculated with these solutions by integrating over all momenta 
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with the correct distribution functions. Thus, the vector and scalar densities of the fermions {i = n, p, t, h) are given 
by 

with degeneracy factor gi = 2 and the energy 



ei{k) = ^Jk^ + K - llif . (59) 

The Fermi-Dirac distribution for the particle {t] = 1) and antiparticle {r] = —1) contributions is defined by 

fnk) = {exp[rJPiEV-^,,)] + l}-' (60) 

where (3 = 1/T and E]'{k) = T,io + riei{k) is the quasiparticle energy. The (relativistic) chemical potential of a particle 
i is denoted by /Xj. The densities of the bosons (i = d, a) are obtained from 

/(fik 
j^^biik)+hi (61) 

„ f k Tfij S7',,._„ , , 



with the Bose-Einstein distribution 



bi{k) = {e^p[(}{E+-i^i)]-l} ' (63) 



and degeneracy factors rjd = 3 and = 1, respectively. A possible contribution to the densities from particles that 
are condensed in the ground state is denoted by hi and n?. In homogeneous and isotropic matter these two are 
actually identical. For a system of nucleons and clusters in chemical equilibrium, the (relativistic) chemical potential 
of a cluster i is determined by 

Hi = NiHn + ZiHp . (64) 

Thus, there are only two independent chemical potentials. 

For given total baryon number density n, asymmetry S and temperature T, the coupled field equations of the 
generalized RMF model are solved selfconsistently. This procedure yields the chemical potentials of neutrons and 
protons that determine the densities of all particles. Finally, all thermodynamical quantities, that are specified in the 
following subsection, can be calculated. 

B. Thermodynamical quantities 

The energy density e and the pressure p are derived from the energy-momentum tensor T'^'' with the results 



12 9iJ fieiik) +Y^9ij hei{k) + hi {rm - E^) 



i=n,p,t,h v i=d,a 



(65) 



and 



3 T. iT-'n -\ (t)3 + 3 .5 (1)3 ^'ei{k) 

+ [rin + rip) [r^wo^w + '^'pPonp - F'^an^] - ^ [m^cr^ - to^Wq - ™p 



(66) 
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with po = (^0)3- The condensed bosons do not contribute to the pressure but to the energy density. The entropy 
density s can be extracted from the grand-canonical potential density ijj(r, /i„, /ij,) = —p as s = — {duj/dT)\^ . 
However, it is more practical to use the thermodynamic relation 

e = Ts - p + ^ fiiUi (67) 

i—n,p,d.t,h,a 

corresponding to the Hugenholtz-van-Hove theorem. After partial integration the standard result 

« = - E ffJ7SjE[/'l^/' + (l-/')l^(l-/')] (68) 

i=n,p.,t,h y ' n 

- E f W TT^I [h\nb,-{l + h)\n{l + b,)] 

i=d,a •' ^ 

is obtained. The thermodynamical pressure 

(69) 



d_n 

dn \ n 



T,S 



calculated from the free energy density / = £ — Ts is identical to the pressure ((66)) in the fieldtheoretical approach. 
Since the energy density in the RMF model contains the contribution of the rest mass of the particles, it is convenient 
to define the internal energy per nucleon as 

EA{n, S,T) = ^ [sin, S, T) - n\°^mn - ^'m^] (70) 

and correspondingly the free energy per nucleon 

FAin,S,T)^EA-T-, (71) 
n 

where the rest mass has been subtracted. We emphasize that both the generalized RMF model and the QS approach 
are thermodynamically consistent. 



C. Dissolution of clusters 



In the generalized RMF approach, the quasiparticle energy shift of the clusters contains the self-energy shift due 
to the mean fields of the mesons. The effective-mass shift is already included in the relativistic approach and the 
Coulomb shift is neglected. For the binding energy shift ABi only the effect of the Pauli blocking is considered. The 
dependence of the Pauli shift on the cluster momentum relative to the medium is neglected since the introduction of a 
momentum dependence in the RMF Lagrangian is non-trivial [46j . With increasing density of the medium, the linear 
approximation of the Pauli shift (j25p in the densities is not sufficient and higher-order terms have to be considered. 
In the present RMF calculation, an empirical quadratic form 



1 



2nO(T) 



SB,{T) 



(72) 



is used, where the abbreviation 



(73) 



and the density scale 



B' 



SBiiT) 



(74) 



for the dissolution of the cluster i with the vacuum binding energy B^ are introduced. The quantity SBi{T) is given 
by (5i?f (T, 0) in eq. for deuterons and in eq. ^7} for the other clusters. In the limit T ^ the shifts and their 
derivatives remain finite for all clusters. 
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FIG. 1: (Color online) Change of the binding energy Bi — Bf + ABi of the clusters i — d,t,h,a at rest in symmetric nuclear 
matter due to the binding energy shift ABi as used in the generalized RMF model as a function of the total nucleon density 
n = nJi°* + 71*°* of the medium for various temperatures T. 



The total binding energy of a cluster i is the sum Bi = + ABi of the experimental binding energy Bf in the 
vacuum [g^I and the binding energy shift ABi that in general depends on the cm. momentum K (see 111 Cp . The 
dependence of the cluster binding energies on the total nucleon density n = + n^f' of the medium is depicted 
for symmetric nuclear matter in Fig. [T] for various temperatures T and clusters at rest. For n — > the experimental 
binding energy is recovered. The density where a cluster becomes unbound, i.e. Bi — 0, increases with increasing 
temperature. This behavior is expected since the Pauli blocking of states is less effective at higher temperatures. Note 
that both free nucleons and nucleons bound in clusters are relevant for the Pauli principle occupying phase space in 
momentum representation. In principle, this is not described by a Fermi distribution but by the bound state wave 
functions. A more exact theory taking this into account is given by the cluster mean-field approximation [s^ . which, 
however, is very complex and has not been solved in general so far. Here we use the approximation of an effective 
Fermi distribution with an effective chemical potential which includes both free and bound state nucleons, as discussed 
following Eq. 1171 It is clearly seen that the weakly bound deuteron dissolves in the medium at much lower densities 
than the more tightly bound a-particle. 

The quadratic form ([72)1 of the binding energy shift ABi predicts a transition of the cluster bound state to the 
continuum at a transition density of n\[T) = (^/S — l) n'^{T) where the cluster binding energy Bi becomes zero, i.e. 
there are only scattering correlations remaining and the energy corresponds to that of a resonance. In case of the 
triton, helion and a-particlc, the resonance energy rapidly moves to larger energies in the continuum with increasing 
density leading to a strong suppression of the cluster fractions. In contrast, the deutcron-like resonance stays closer 
to the threshold causing a much weaker suppression of two-particle correlations at high densities. The change of 
the binding energy from positive to negative values allows to decribe a continuous suppression of the cluster fraction 
with increasing density. A simple neglection of the cluster contribution to the EoS as soon the energy crosses zero 
would lead to an unphysical jump in the cluster density and the thermodynamical properties. In the QS model, the 
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TABLE 11: Masses of the nucleons and mesons in the relativistic mean-field modeL 
particle i neutron proton lo meson a meson p meson 
'rhi [MeV] 939.56536 938.27203 783 546.212459 763 



TABLE III: Parameters of the couplings in the relativistic mean-field model. 



meson i 


ri(nsat) 


ai 




Ci 


di 


uo 


13.342362 


1.369718 


0.496475 


0.817753 


0.638452 


a 


10.686681 


1.357630 


0.634442 


1.005358 


0.575810 


P 


3.626940 


0.518903 









bound state contributions are rapidly cancelled by the continuum contributions, Eqs. (|30l3ip . leading to a more rapid 
suppression also of the deuteron correlations. The differences in the continuum correlations will be seen to have large 
effects in the comparison of the results from the RMF and QS models. 

D. Model parameters 

The generalized RMF model contains several parameters: the masses of the particles, the couplings and binding 
energy shifts with their specific functional dependence on densities and temperature. In the present approach, exper- 
imental neutron and proton masses m„ and m,p are used instead of an average nucleon mass Wnuc — i'm-n + "^p)/2. 
With the experimental binding energies B° in the vacuum from (60i] the masses ([36]) of the clusters are also fixed. For 
the masses of the uj and p meson standard values of previous RMF models are assumed. The mass of the a meson is 
determined from a fit of the RMF parameters to properties of finite nuclei (see below) . The numerical values of the 
nucleon and meson masses are given in table [Til 

The functional dependence of the couplings on the density is described by 

with X = n/usat where a rational function 

l + bi{x + di)'^ 

is used for the isoscalar mesons i — uj,a and an exponential function 

fi{x) = exp[-ai(a; - 1)] 

for the isovector meson i = p. In order to reduce the number of independent parameters the conditions /i(l) = 1 
and /j"(0) = are imposed on the rational function. These functions were introduced in [i^ and are now widely 
used in RMF models with density-dependent couplings. The saturation density rigat, the mass of the a meson mo-, 
the couplings Ti{nsat) and the coefficients a^, bi, Ci and di are found by fitting the properties of finite nuclei (binding 
energies, spin-orbit splittings, charge and diffraction radii, surface thicknesses and the neutron skin thickness of ^'^^Pb) 
in the same way as for the parametrization DD in Ref. [46j . In total, there are ten independent parameters in the fit. 
Numerical values of the coupling parameters can be found in table IIIIl This new parametrization is called DD2 since 
it is a modification of the set DD where the only difference is the use of experimental nucleon masses. 

On the basis of this fit, the saturation density of symmetric nuclear matter at zero temperature is obtained as 
"■sat — 0.149065 fm~'^ with a binding energy per nucleon of —16.02 MeV. The incompressibility turns out to be 
= 242.7 MeV with a derivative K' = -529.8 MeV. See [H,!!!] for the definition of these quantities. These values 
are very reasonable and close to the results of other modern RMF parametrizations. The large negative value of K' is 
a result of the fit to the surface properties of nuclei and leads a rather stiff EoS for symmetric nuclear matter at high 
densities. The small effective Dirac mass at saturation of 0.5625 Wnuc is required in order to get a good description 
of the spin-orbit splittings. This corresponds to an effective Landau mass of m* = 0.6255 mnuc- 

The value and the density dependence of the symmetry energy near the nuclear saturation density rigat are usually 
characterized by the quantities J — i?sym('^sat, 0, 0) and the slope parameter L = 3 dEsym/ dn\^^^ ^ j,^q. With the 
parametrization of the present RMF model, the values J = 32.73 MeV and L = 57.94 MeV are found. The obtained 
symmetry energy at saturation J is fully consistent with all modern RMF parametrizations and expectations. The 
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FIG. 2: (Color online) Proton fraction Xp in symmetric nuclear matter as a function of the total density n for various 
temperatures T in the generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The result 
of the NSE calculation with light clusters is denoted by thin solid lines for low densities. 



rather small slope coefficient L is a consequence of fitting the neutron skin thickness of ^"^Pb that is not known 
precisely so far. Similar low values for L are found in other contemporary RMF parametrizations with density 
dependent couplings or with extended non-linear meson self-interactions. Older non-linear RMF models were not able 
to give a reasonable value of the neutron skin thickness with values for L in excess of 100 MeV. They displayed a much 
stiffer symmetry energy because of a restricted form for the isospin dependence of the interaction. The symmetry 
energy of nucler matter parameters at subsaturation densities is one of the central results of this work and will be 
discussed in section IVll 



IV. PROPERTIES OF SYMMETRIC NUCLEAR MATTER WITH LIGHT CLUSTERS 

The appearance of light clusters in nuclear matter at densities below saturation affects the composition and the 
thermodynamical properties of the system. In this section we will compare the results of the quantum statistical 
approach (QS) with the generalized relativistic mean-field model (RMF) in reference to the nuclear statistical equi- 
librium (NSE) model, which gives the correct behavior in the limit of small densities. We will start to discuss the 
composition of the system, which shows most directly the differences of the models. In all figures in this section we 
consider isothermes in symmetric nuclear matter as a function of the total baryon density for temperatures between 
2 MeV and 20 MeV in steps of 2 MeV keeping the same color code. This representation immediately allows to study 
the systematic evolution of the various properties. 



A. Composition 

We start to discuss the composition of the system, i.e. the fractions Xi = AiTii/n of the various particle species i 
as calculated in the RMF and the QS models. We first describe the results and then attempt to give an explanation 
of the differences between the models at the end of the subsection. 

In Fig. [2] the density dependence of the free proton fraction, i.e. of protons not bound in a cluster, is shown for 
different temperatures in symmetric nuclear matter (which in this case is nearly identical to the neutron fraction). 
The two approaches are compared to the NSE model up to densities n < 2 - 10~^ fm~'^. For finite temperatures T the 
free proton fraction in symmetric nuclear matter always approaches the value 0.5 for n ^ fm^"^. It first decreases 
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FIG. 3: (Color online) Cluster fractions Xi in symmetric nuclear matter as a function of the total density n for various 
temperatures T in the generalized RMF model (thick solid lines). The result of the NSE calculation with light clusters is 
denoted by thin solid lines for low densities. See Fig. |2] for the color code. 



with increasing densities because of the formation of clusters, but then increases again because the clusters dissolve 
at higher densities. Eventually the system becomes homogeneous again, and the nucleon fractions attain the value 
0.5. For densities below n ^ 10"'* fm~'^ the fraction of free protons in both models is very well described by the NSE 
result since here mean-field effects and changes of the cluster properties are practically negligible. With increasing 
density, the NSE proton fraction approaches zero asymptotically irrespective of the temperature, i.e. all protons are 
predicted to be bound in clusters. This unphysical result does not occur in both the RMF and QS approaches. 
Instead, the clusters dissolve at high densities leading to free protons and neutrons at high densities, i.e. the correct 
limit is obtained. 

The two approaches, RMF and QS, generally show a similar behavior. However, in the transition region where 
the clusters dissolve, there are significant differences, which are more pronounced at high temperatures. In the QS 
approach the minimum of the free proton fraction increases with the temperature monotonously, and the minimum 
position is slightly shifted to higher values in density but stays close to 10"^ fm"^. Essentially all protons are free 
at saturation density independent of temperature. In the RMF model the behavior at temperatures above w 8 MeV 
is different. The minimum of Xp starts to decrease again and the minimum position moves considerably to higher 
densities with increasing temperature, such that the model predicts that some protons are still bound in clusters even 
at saturation density. 

The information on the fractions of the different light clusters, deuteron to *He, is shown in Figs. [3] for the RMF 
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FIG. 4: (Color online) Cluster fractions Xi in symmetric nuclear matter as a function of the total density n for various 
temperatures T in the quantum statistical approach (thick solid lines). The result of the NSE calculation with light clusters is 
denoted by thin solid lines for low densities. See Fig. |2] for the color code. 



model and in Fig. |3]for the QS model, again in comparison with NSE. At zero temperature and density the cluster 
fractions are simply determined by the binding energies. Thus in symmetric nuclear matter all nucleons are bound 
in alpha particles. With increasing temperature the cluster fractions show the complementary behavior compared to 
the free nucleon fractions, discussed above. At low densities and higher temperatures all cluster fractions are small 
decreasing with temperature but increasing with density. For a fixed temperature first the deuteron appears, then 
the three-body bound states and finally the a particle. The distribution between the clusters changes with increasing 
density of the medium such that correlations between more and more nucleons become important. As a consequence, 
the fractions of the lighter clusters decrease whereas the a-particle fraction still increases. In the QS model (but not in 
the RMF model) at higher densities the more strongly bound a-particle are formed more frequently than deuterons. 
The same trend is observed for the NSE calculation. However, for densities beyond nuclear saturation, the NSE 
predicts that all nucleons in symmetric nuclear matter would be bound in clusters, which again is unphysical since in 
this model the medium modification and the eventual disappearance of the clusters is not taken into account. 

Again there are significant differences in the two approaches. In the RMF model the cluster fractions along 
isothermes rise and fall monotonously with a single maximum. In general they stay below the fractions predicted 
by the NSE calculation, since the clusters are less bound inside the medium. In the QS approach, on the other 
hand, at lower temperatures, the fractions of deuterons, tritons and helions exhibit a sudden drop around densities of 
10~^ fm~'^, which is accompanied by an increase in the a-particle fraction as compared to the NSE result. Generally, 
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the range of densities where the clusters disappear is more confined in the QS relative to the RMF model. 

Many of the differences in the behavior of the models can be traced back to the different treatment of the deuteron 
correlations. The QS approach takes the continuum contributions explicitly into account, which effectively reduces 
the strength of the two-body correlations, as is, e.g., seen in Eqs. (|30p . PT|) . In RMF, on the other hand, the 
deuteron correlations are represented by a single state that slowly moves to higher energies in the continuum, and 
thus two-body correlations are overestimated. Also in NSE the deuteron correlations are overestimated, since there 
is no continuum contribution. This has an strong effect on the deuteron fraction and - in competition - also on the 
fractions of nucleons and other clusters. We see that in the RMF model the deuteron fractions at higher densities are 
generally larger than in the QS model, as seen by comparing the upper left panels of Figs. [3] and[4l This, in turn, 
has a strong effect on the a fractions, which are much lower in the RMF model as seen in the lower right panels of 
these figures, but also leads to the decrease of the free nucleon fraction at higher densities, as seen in Fig. [51 These 
effects can also be seen in the comparison with the NSE limit, where they lead to larger deuteron fractions relative 
to the QS model. Only at temperatures smaller than the deuteron binding energy, the bound state is the dominating 
two-body correlation. A similar effect occurs in the virial description of matter at low densities that is encoded in the 
temperature denpendence of the second virial coefficient. For heavier clusters the influence of the continuum on the 
fractions is much less pronounced due to their larger binding energies. 

As we remarked above, the QS model shows a particularly enhanced alpha particle fraction at the higher densities, 
as seen in the increase above the NSE limit at densities around n ~ 10~^ fm~'^ in Fig. 01 This has much of the 
appearance of an onset of an alpha particle condensation. As a consequence, the fractions of the other clusters show 
a dip around this density, and there are also consequences in the thermodynamical quantities, as seen below. All 
these effects are not present in the RMF model. However, mean field contributions from the rather substantial cluster 
fractions in this density range are not taken into account in the QS model, the effect of which needs to be further 
investigated. It is important to note that the a particle is usually not the most frequent cluster and that there are 
substantial contributions from the deuteron, the triton and the helion at intermediate densities and temperatures 
that are not considered in the EoS of Lattimer and Swesty [l^ or Shen, Toki et al. [2§|. Furthermore, the excluded- 
volume mechanism to suppress the formation of clusters at high densities does not take into account any temperature 
dependence in this process which is clearly present in our more microscopic models. 

B. Thermodynamical quantities 

In this subsection we discuss the thermodynamical quantities for symmetric nuclear matter. The essential effects of 
the formation of clusters in the various models are given by the pressure p as a function of the total baryon number 
density density n. For a better representation, we depict in Fig. O the ratio p/n as a function of n. The left and 
right panels of the figure show the results of the RMF and the QS models, respectively, with thick lines. The thin 
lines in both panels represent a NSE calculation with neutrons, protons and the light clusters d, t, h and a. They are 
shown only for densities below 2 • 10~^ fm~'^ because the contribution of heavier clusters can be substantial at higher 
densities, at least for low temperatures, and the NSE becomes unrealistic. 

At very low densities, all models approach the ideal gas limit with p/n = T since the cluster fraction is very small 
and nuclear matter is composed primarily of neutrons and protons (see Subsect. IIV A[) . Relativistic and mean field 
effects are not important here. With increasing density the NSE calculation exhibits a reduction of the pressure 
relative to the ideal gas that is caused by the formation of light clusters. Both the RMF and the QS models follow 
this trend but there are considerable differences. The QS results stays closer to the NSE calculation whereas RMF 
shows a deviation already at lower densities around 10^"^ fm. However, both models predict a stronger decrease of the 
pressure than the NSE, because at these densities the changes of the cluster properties and the mean-field effects are 
already effective. An exception is seen for the RMF approach at temperatures below approx. 6 MeV where it shows 
an increase relative to the NSE. 

In the density range from 10"^ to 2 • 10~^ fm~^ the ratio p/n displays considerable structure and passes through 
at least one minimum, eventually with negative pressure, before it rises sharply for densities above nuclear saturation 
{p/n = 0) and the matter becomes very incompressible. In this region, we observed the most pronounced differences 
in the behavior of the RMF and QS models with respect to the composition in Subsect. IIV Al They are generated 
partly by differences in the description of the cluster correlations and partly by a different treatment of the mean-field 
effects. In the QS approach, the strengths of the mean fields and thus the quasi-particle energies and mass shifts are 
independent of the composition of the system since they are taken in parametrized form from the RMF calculation 
without clusters. In contrast, in the RMF approach the additional contributions in the source terms of the meson fields, 
see Eqs. (H5|) and ([^ . depend on the cluster densities. They describe the back reaction of the cluster formation on the 
medium (apart from the contribution of the bound nucleons) and are inevitable for the thermodynamical consistency 
of this model. In the density region where the cluster fraction is substantial, the additional terms lead to a sizable 
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FIG. 5: (Color online) Ratio pressure over density p/n of symmetric nuclear matter as a function of the total density n for 
various temperatures T in the generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The 
result of the NSE calculation with light clusters is denoted by thin solid lines for low densities. See Fig. [2] for the color code. 



modification of the vector meson fields. In the density range around n ^ 10^^ fm^"^ the pressure is lower in the QS 
relative to the RMF model. This correlates with increased alpha formation (a "condensation" ) in this model relative 
to RMF, as was seen in Fig. H) 

The relativistic baryon chemical potentials /i = (/j,p + /i„)/2 of the RMF and the QS approaches are displayed in Fig. 
[6] in comparision with the NSE calculation. They refiect the behavior observed for the pressure of the various models. 
At very low densities, the three calculations agree with each other and the ideal gas dependence of the chemical 
potential /i = m + Tin on temperature and density. With increasing density, the chemical potential of both 

the RMF and QS calculations are below the NSE results, except for the RMF model at low temperatures as already 
seen for the pressure. Here, the chemical potential rises above the NSE result with increasing density and a pronounced 
maximum occurs. In the QS model, the chemical potential generally stays below the NSE result. Deviations from 
the NSE predictions appear already at lower densities in the RMF relative to the QS model. In general, however, the 
differences between two approaches are less obvious in the chemical potential than in the pressure. In both models 
the low-density and the high-density limit are correctly described. 

The density dependence of the free energy per nucleon Fa, Eq. ([7T|) . is shown in Fig. [7] and closely follows the 
density dependence of the chemical potential. In the QS approach. Fa is obtained by a simple direct integration of 
the chemical potential, cf. Eq. (fTU]) . Both models agree perfectly with the NSE result for densities below 10"'' fm~^. 
At the nuclear saturation density, the free binding energy per nucleon approaches the local minumum at —16 MeV in 
the limit T ^ MeV as expected. 

The results for the internal energy per nucleon Ea are depicted in Fig. [51 Here the differences between the RMF 
and QS results are less pronounced than in the chemical potential ^ or the free binding energy per nucleon Fa- In 
contrast to the latter quantities, Ea increases with temperature at a given density. For n ^ the classical limit 
Ea 3T/2 of an ideal gas is approached. Slight deviations stem from the use of relativistic dispersion relations. 

The entropy per nucleon Sa = s/n, shown in Fig. [9l generally decreases monotonously with increasing density n 
except for a small range in density at low temperatures when the light clusters suddenly dissolve in the medium. At 
low densities the entropy per nucleon in the RMFand QS models approaches the NSE result and at very low densities, 
a dependence Sa oc — ln(n) -I- const, is found consistent with the behavior of an ideal gas. The occurence of clusters 
is found to lead to a reduction of the entropy per nucleon as compared to pure neutron-proton matter. 
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FIG. 6: (Color online) Relativistic baryon chemical potential fi = (/Xp + ^in)/2 of symmetric nuclear matter as a function of 
the total density n for various temperatures T in the generalized RMF model (a) and the QS approach (b) with light clusters 
(thick solid lines). The result of the NSE calculation with light clusters is denoted by thin solid lines for low densities. See 
Fig. [2] for the color code. 




FIG. 7: (Color online) Free energy per nucleon Fa of symmetric nuclear matter as a function of the total density n for various 
temperatures T in the generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The result 
of the NSE calculation with light clusters is denoted by thin solid lines for low densities. See Fig. [2] for the color code. 
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FIG. 8: (Color online) Internal energy per nucleon Ea of symmetric nuclear matter as a function of the total density n for 
various temperatures T in the generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The 
result of the NSE calculation with light clusters is denoted by thin solid lines for low densities. See Fig. [2] for the color code. 




FIG. 9: (Color online) Entropy per nucleon Sa of symmetric nuclear matter as a function of the total density n for various 
temperatures T in the generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The result 
of the NSE calculation with light clusters is denoted by thin solid lines for low densities. See Fig. [2] for the color code. 
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FIG. 10: (Color online) Pressure (a) and relativistic baryon chemical potential (b) of symmetric nuclear matter as a function of 
the baryon number density in the generalized RMF model with clusters (solid lines) and without clusters (dashed line) taking 
the phase transition into account. See Fig. [5] for the color code. 

V. LIQUID-GAS PHASE TRANSITION 

In a system of given total neutron and protons number densities, n^°* and rip"*, and temperature T, the corresponding 
thermodynamical potential, i.e. the free energy density np°',T), should be minimized. However, this procedure 

does not necessarily give the correct equilibrium state. Thermodynamic laws require that the free energy density 
is a convex function in the variables n^°*, n^"* and T to assure the stability of the system. This condition leads 
to the occurence of phase transitions. The coexistence region of different phases in thermodynamical equilibrium is 
separated from the region of a single phase by the so-called binodal surface. It can be obtained for a given temperature 
by a general Gibbs construction where the values of the intensive variables pressure p and the chemical potentials 
of protons /ip and neutrons /i„ in the two phases have to be identical, see, e.g. [6ll . m^ . Thus, the binodal surface 
is found from a global criterion in contrast to the spinodal surface that defines the boundary of local instability of 
the system, which occurs, e.g., when the compressibilities of the system become negative. In general, the spinodal 
is enclosed by the binodal and both become identical along critical lines. In the case of symmetric nuclear matter, 
the problem becomes one-dimensional and a usual Maxwell construction for the phase transition is sufficient with 
constant pressure p and baryon chemical potential /i for densities inside the coexistence region. In the general case 
of asymmetric matter, the pressure and chemical potentials do not stay constant. In addition, the low-density phase 
has a larger isospin asymmetry than the coexisting high-density phase. 

Nuclear matter as considered here is an idealized physical system where the Coulomb interaction is neglected and 
charge neutrality is not demanded. The phase transition boundary as constructed by the above standard procedure 
only gives a first indication of where the formation of inhomogeneities and heavy nuclei occurs. In more realistic 
calculations with Coulomb interaction the system has to be globally charge neutral and at least the contribution of 
electrons has to be considered, e.g. in applications of the EoS to astrophysics. We will not follow up on this issue in 
the present work. For references to work in this direction, see our discussion of pasta phases in the Introduction. 

The results for the pressure p, Fig. [SI and the baryon chemical potential fj,, Fig. [SI allow to construct the phase 
transition in symmetric nuclear matter from the gas phase with nucleons and light clusters at small densities to the 
liquid phase of pure nucleon matter at high densities. In the following we consider the example of the RMF model. 
In Fig. [TUj the pressure p and the chemical potential fi with phase transition are shown for the cases without and 
with light clusters. Inside the coexistence region of the two phases pressure and chemical potential stay constant 
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FIG. 11: (Color online) Binodal line enclosing the liquid-gas coexistence region in the pressure-density diagram (a) and phase 
transition line in the temperature-chemical potential diagram (b) for symmetric nuclear matter in the RMF model without 
light clusters (blue lines), with light clusters (red lines) and in the QS approach (green lines). The corresponding critical points 
are denoted by circles. 



as characteristic for a Maxwell construction. As the temperature approaches zero, the density range of this region 
becomes larger with the limits p — > MeV fm^^ and /j, ^ m — 16 MeV. The differences between the calculations 
without and with clusters increase at higher temperatures with a larger coexistence pressure in the latter case. Finally, 
the critical temperature is reached beyond which no phase transition occurs any more. 

At low temperatures the binodal line at the lower boundary of the coexistence region is reached already at very 
small densities. At this point the cluster fraction is still rather small, and, consequently, the appearance of clusters 
has little effect on the determination of the phase boundary for low temperatures. Larger effects of the clusters on the 
phase transition are observed only at higher temperatures where the cluster fraction at the lower density boundary 
reaches larger values. The binodal on the high density side of the phase coexistence region is hardly affected when 
clusters are considered in the calculation because here the matter is essentially composed of free nucleons. 

Results as in Fig. [10] can be used to extract the binodal line enclosing the phase coexistence region in the p — n 
diagram and the phase transition line in the T — ^ diagram. These are shown in the left and right parts of Fig. 1111 
respectively, for the RMF model without and with clusters and the QS approach. The occurence of light clusters 
in the system narrows the width of the coexistence region and shifts the maximum to a higher density and larger 
pressure in the generalized RMF model as compared to the calculation without clusters. In the RMF model the 
critical temperature Tc increases from 13.72 MeV without clusters to 15.12 MeV with clusters. The corresponding 
values of the critical density Uc are 0.0452 fm^'^ and 0.1018 fm~'^, and of the critical pressure Pc 0.1781 MeV fm~'^ 
and 0.9029 MeV fm~^, respectively. Thus, in the RMF model the position of the critical point shifts rather drastically 
when light clusters are considered. There is also a marked effect on the phase transition line in the T — fj, diagram. 
In RMF without clusters, the chemical potential decreases monotonously with increasing temperature until it reaches 
the critial point at a critical chemical potential of = 914.48 MeV. At low temperatures the phase transition line 
with clusters follows closely the line without clusters. However, towards the critical point the transition line with 
light clusters has as S shape. At temperatures T above « 8 MeV the curve finally bends to higher chemical potentials 
ending at the critical point with = 918.9 MeV. 

In the QS approach a very different trend is found. The binodal line exhibits a characteristic dip on the low-density 
side that is related to the enhancement of the a-particle fraction, cf. Fig. ID The critical point moves to a smaller 
temperature of Tc = 12.1 MeV and a slightly higher chemical potential of /ic — 915.61 MeV as compared to the 
RMF calculation without clusters. The critical pressure hardly changes. The observed differences between the RMF 
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FIG. 12: (Color online) Free energy per nucleon Fa (left) and internal energy per nucleon Ea (right) as a function of the 
asymmetry 5 for constant total density n = 0.001 fm~^ for various temperatures T without light clusters (dashed lines) and 
with light clusters (solid lines). See Fig [2] for the color code. 



and the QS approaches are again related to the fact that many-body correlations at high temperatures, especially of 
the deuteron, are overestimated in the former. The behavior of the RMF approach must be considered unphysical 
in this respect, because one would expect, that correlations decrease the chemical potential and the pressure of the 
phase transition. In an improved version of the generalized RMF approach which takes the continuum contributions 
explicitly into account, it is expected that the critical point moves to lower temperatures. 



VI. SYMMETRY ENERGY 



The interanl energy per nucleon Ea of asymmetric nuclear matter, as defined in Eq. (j70p . can be expanded for 
given density n and temperature T in powers of the asymmetry S 

EA{n,5,T) = EA{n,0,T) + Esy^{n,T)S^ + . . . . (78) 

The symmetry energy is the coefficient of the first term in the expansion that depends on the asymmetry S quadrati- 
cally. Thus, it is defined as the second derivative 

(79) 

<5=0 

Usually, the dependence of EA{n, S, T) for the complete range of asymmetries is quite well approximated by a quadratic 
function. In this case the symmetry energy also represents the difference between the binding energy per nucleon of 
neutron matter and of symmetric nuclear matter. At finite temperatures one has to distinguish between the internal 
symmetry energy Esymin,T) and the free symmetry energy Fsymiji^T) that is similarly defined. 

However, the quadratic approximation is not valid in general, especially for a system with cluster correlations at 
densities n below the nuclear saturation density. This is clearly seen for our models in the dependence of Ea and 
Fa on 6 for given n and T. In Fig. [T^ the free binding energy per nucleon (left) and the internal binding energy 
per nucleon (right) are shown as a function of the asymmetry for a fixed density n = 0.001 fm~^ and different 
temperatures. Without clusters, the dependence on 5 is rather weak and can be well decribed by a parabolic function 
irrespective of the temperature. With clusters, the binding energies per nucleon are substantially lowered around 
(5 = 0, i.e. symmetric nuclear matter, particularly for low temperatures. The system gains additional binding energy 
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FIG. 13: (Color online) Internal symmetry energy i?sym as a function of the total density n for various temperatures T in the 
generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The result of the RMF calculation 
without clusters is denoted by thin solid lines for low densities. See Fig [5] for the color code. 

by forming clusters. Large deviations from a global quadratic dependence develop and especially at low temperatures 
the parabola changes to a triangular shape, most obviously seen for the internal energy per nucleon. This behavior 
leads to a very large and thus not very meaningful symmetry energy in the limit T — > MeV when the conventional 
definition (j79p is used. Hence, a more appropriate characterization for the symmetry energy is required. A reasonable 
choice is given by the finite difference formula 

-BsymKT) = ]^[EA[n,l,T) ~2EaMT) + EA{n,-l,T)] (80) 

that is identical to (|79|) for an exact quadratic dependence of the internal binding energy per nucleon on 6. This 
modified definition gives a good measure of the binding energy differences between neutron, proton, and symmetric 
nuclear matter. A corresponding equation defines the free symmetry energy Fsym{'n,T). 

The density dependence of the internal symmetry energy i^sym and of the free symmetry energy i^sym as defined 
above is presented in Figs. [T51 and [HI respectively, as a function of density for various temperatures. In these figures, 
the results of the RMF and QS approaches with clusters are compared to the RMF result without clusters. Just as for 
the free energy. Fig. [7l and the internal energy. Fig. [8l of symmetric nuclear matter, the free and internal symmetry 
energies exhibit a different behavior for n — > 0. In all models, i?sym approaches zero in this limit, but i^sym converges 
to Tin 2. 

In the RMF model without clusters both internal and free symmetry energies rise continuously with increasing 
density, and i?sym is almost independent of T. When the formation of clusters is taken into account, the internal 
symmetry energy increases substantially at low densities, see Fig.[T31 This behavior is caused by the additional binding 
of symmetric nuclear matter, that was already seen in Fig. 1121 which is particularly pronounced at low temperatures 
with a large cluster fraction. The density dependence of the internal symmetry energy is rather different for the for 
the RMF and QS approaches in a region near 10~^ fm^'^. As was discussed in Subsect. IIV Al in the RMF calculation 
the fraction of deuterons is enhanced in the transition region whereas three- and four-body correlations are suppressed 
as compared to the QS approach. Correspondingly, symmetric nuclear matter is less bound here and the symmetry 
energy is reduced relative to the QS approach. Since two-body correlations survive in the RMF model to too high 
densities for larger temperatures, the internal symmetry energy is also lowered considerably when the total baryon 
number density approaches the saturation density. In the QS approach, Ea is much closer to the RMF result without 
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FIG. 14: (Color online) Free symmetry energy Faym as a function of the total density n for various temperatures T in the 
generalized RMF model (a) and the QS approach (b) with light clusters (thick solid lines). The result of the RMF calculation 
without clusters is denoted by thin solid lines for low densities. See Fig. [5] for the color code. 



clusters with increasing density. The free symmetry energy, see Fig. [HI shows essentially the same features as the 
internal symmetry. However, they are less pronounced because of the different low-temperature limit. 



VII. COMPARISON WITH OTHER APPROACHES 



We discuss light cluster abundances obtained from our approaches in comparison with results from other models, 
which appear in the literature and which have been discussed in this paper. These are: the EoS of Shen et al. [2^, an 
EoS based on a virial expansion as presented in [s^ . and a nuclear statistical equilibrium calculation that takes into 
account the ground states of all nuclei in the atomic mass evaluation (AME 2003) [60|]. Fig. [T5l shows the a-particle 
fraction in symmetric nuclear matter at four different temperatures as a function of density. We have not explicitly 
included the results of the model of Lattimcr and Swesty 28] in this comparison. They use the same excluded volume 
prescription as Shen et al. [29] to take into account the medium dependence of the a-particle fraction. Therefore, the 
behavior with density of their results is similar to that of Shen et al. , even though absolutely there are differences due 
to a different treatment of the heavy particle fraction. A direct comparison of of Shen et al. and Lattimer et al. 



[28l |. can be found in Ref. [31 1 



At low densities all models show a decreasing a-particle fraction with decreasing total baryon number density. How- 
ever, there are some subtle differences. The virial expansion (black dash-dotted lines) takes scattering contributions 
to the second virial coefficient into account. These are important at higher temperatures and cause a slightly higher 
prediction for Xa at small n. The effect of cluster dissociation is absent in this model and results in an monotonic 
increase of the cluster fraction with increasing density, which is unphysical since one does not expect to see clusters 
in nuclear matter at saturation density. 

The NSE approach (green dotted lines) compared to the virial approach shows a reduction of the a-particle fraction 
at all densities, due to the formation of heavier clusters, a-particles survive even beyond the nuclear saturation 
density, since medium modification and dissolution of clusters are absent in this model, which also fails to describe 
the transition from clusterized matter to cluster-free nuclear matter at high densities. 

The Shen et al. EoS (blue dashed lines) neglects the contribution of deuteron, triton and hclion clusters. This 
leads to an overestimation of the a-particle fraction at low densities where actually lighter clusters dominate the 
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FIG. 15; (Color online) Comparison of Q-particle fractions in symmetric nuclear matter as a function of the density at four 
temperatures for the virial expansion (black dashed-dotted lines), NSE (green dotted lines), the EoS of Shen et al. 29] (blue 
dashed lines), the generalized RMF model (red solid lines) and the QS approach (orange dashed lines). Note the different scales 
on the a;- axes. 



composition. In contrast to the above two approaches, there is a steep decrease of the a-cluster abundance when 
approaching the saturation density. However, one notes some irregularities which are understood as an effect of the 
excluded volume approach when the closest packing is reached. 

The generalized RMF model (red solid lines) developed in this work describes the decrease of the a-particle fraction 
at high densities by a reduction of their binding energy due to the Pauli blocking which leads to the Mott effect 
for vanishing binding. The maximum cluster density is reached around the Mott density. Due to the presence of 
strong correlations in the scattering state continuum which are effectively represented by one resonance, there is a 
nonvanishing cluster fraction above the Mott density. Among all models presented in the comparison of Fig. 1151 the 
generalized RMF approach shows the strongest reduction of the a-cluster fraction. 

In the QS approach the behavior of the density dependence of Xa is similar to that of the generalized RMF model, 
but has for all densities a higher a-cluster fraction, which is accompanied by smaller light cluster abundances. As 
discussed in Subsect. IIV Al the difference in Xa between the two models is mainly due to the overemphasis of the 
deuteron correlations in RMF, which suppresses the a-particle fraction. On the other hand, the RMF takes into 
account the back-reaction of the cluster formation on the mean-field which is missing in the QS approach. 
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The difference between tfie two approacfies of this paper, shown in Fig. [15] in direct comparison, also gives an 
indication of the effects of possible improvements in the models, apart, of course, from the inclusion of heavier 
clusters, resp. of nuclei embedded in a clusterized gas. On the other hand, the synopsis of Fig. [15] demonstrates the 
advantage of systematic many-body approaches to the description of cluster formation over alternative approaches 
which lack a microphysical mechanism to account for the cluster breakup, but which are extensively used, e.g., in 
astrophysical applications. 

VIII. CONCLUSIONS 

Up to now there exist different strategies to model the EoS of nuclear matter and, in particular, to extract the 
symmetry energy: (1) phenomenological density functional methods, such as non-relativistic Skyrme or relativistic 
RMF functionals, (2) effective-field theory approaches based on density functional or chiral perturbation theory, or 
(3) ab initio approaches, such as Brueckner-type methods, variational calculations or Green function methods. A 
recent overview over these methods with references can be found, e.g., in Ref. [52l |. These methods are developed 
to obtain reliable equations of state for nuclear matter for a range of densities and asymmetries. In particular, the 
density-dependent RMF model can be considered phenomenologically a very useful approach in this respect since it 
has been applied with great success in the simultaneous description of cold nuclear matter in finite nuclei and compact 
stars as well as of hot nuclear matter in heavy-ion collisions and supernova explosions. However, all these approaches 
fail in the low-density limit, where cluster formation becomes essential. In this region, simple approaches which take 
clusters into account, such as the NSE or the Beth-Uhlenbeck formula of the virial expansion, show that clusters give 
a substantial contribution to the composition and the thermodynamic properties. On the other hand, models like the 
NSE or the virial expansion fail at higher densities where in-medium effects become important leading to a dissolution 
of clusters and the transition to cluster-free nuclear matter. 

Here we propose, for the first time, a unified treatment that takes both limits into account describing the smooth 
transition from clusterized matter at low densities to pure nucleonic matter at high densities. We thus suggest a 
symbiotic framework which combines the merits of a QS approach in describing cluster properties in a medium with 
those of the RMF approach to model nucleonic self-energy effects, resulting in two hybrid approaches to the problem 
of cluster formation and dissolution. On the one hand we use a microscopic QS approach to describe the medium 
modification of cluster binding energies due to Pauli blocking. In this approach we employ nucleon self-energies taken 
in parametrized form from a very recent version of the density-dependent RMF model. On the other hand, the results 
of the QS approach for the cluster properties inside the medium are incorporated into an effective hadronic field 
theory using a Lagrangian formulation. Such a theory already has the correct high-density behavior as deduced from 
comparisons to heavy- ion collision experiments. We improve the low-density behavior by explicitly including light 
clusters, such as deuterons, tritons, helions and a-particles as explicit degrees of freedom taking into account the 
medium modification of their binding energies from the results of the QS approach. 

In the numerical evaluation of both hybrid approaches developed in this work, we find that well-defined clusters 
appear only for densities below approximately 1/10 to 1/100 of the saturation density and get dissolved at higher 
densities. A direct confirmation of the given approach can be obtained from a comparison with recent results from 
heavy- ion collisions at low-energies These investigations indicate larger values for the symmetry energy in com- 
parison with the mean field results at low densities, which seem to be in agreement with our findings ^65j . 

Realistic approaches to the clustering in low-density nuclear matter should include excited states such as resonances 
and also the contribution of the continuum of scattering states. This can be done for the second virial coefficient, as 
demonstrated by the generalized Beth-Uhlcnbeck approach [4l| . In this way one can also reach the exact low-de nsity 
limit of the virial approach [sTj, jH, ^ . Compared with the approaches of Lattimer et al. [1^ and Shen et al. [29| . 
who employ phenomenological concepts such as the excluded volume, we have given in this work a description of 
medium effects on the clusters and of their breakup as a result of the fundamental Pauli principle. 

The extension of the present framework to larger clusters beyond the a-particle is straightforward along the lines 
given in this work. For the evaluation of self-energies and Pauli shifts of A-particle clusters in nuclear matter, see [i^. 
The generalization of the given approach to account for clusters of arbitrary size would lead to an improvement in 
the low-density limit when comparing the nuclear statistical equilibrium as used, e.g., in multifragmentation models 
[63l [g^ I . One can alternatively also introduce the formation of heavier nuclei [i^ in the presence of a nucleon and 
cluster gas, in a similar way as it was done in the Thomas-Fermi approximation in the Shen et al. approach. This 
will be relegated to a subsequent paper. We restricted our present work to that region of the phase diagram where 
heavier clusters with A > A are not relevant. 

In our models, we constructed the phase transition from clusterized, gaseous low-density matter to a cluster-free 
nuclear liquid at high densities. The coexistence region gives a hint about the range in temperature and density where 
the occurrence of inhomogeneities and the formation of heavy clusters become relevant. An issue for future investiga- 
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tions is an improved description near the phase transition taking into account effects of the Coulomb interaction and 
charge screening. 

We are able to give the composition and the thermodynamic quantities in a large region of densities, temperatures 
and asymmetries as they are required, e.g., in supernova simulations. We did not consider contributions from, e.g., 
electrons, neutrinos or photons, to the thermodynamical quantities. In astrophysical applications of the EoS, they 
have to be included. They will modify the properties of the system and affect, in particular, the occurence of 
inhomogeneities and of the liquid-gas phase transition. 

As a long-range objective, we aim at a unified description of nuclear matter from very low density to, eventually, 
the deconfinement phase transition that is based on a more microscopic and self-contained description than previous 
approaches to the EoS which are used up to now in astrophysical models. 
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APPENDIX A: LOW-DENSITY EXPANSION 

A density-dependent RMF model was considered in (46j . The following low-density expansions are derived form 
this model and reproduce the DD-RMF results below the baryon density n < 0.2 fm~'^ within 0.1 %; variables are 
the total baryon density n = + in units of fm~'^, the asymmetry 6 — — + "-^T*) the 

temperature T in MeV. The scalar field in MeV is given by 

E„,p(r, n, S) ^ ns[(4524.13 - 6.926 T) - UM57S^/4 + 0.833943 ^'*/16 - 9.00693 6^64] 
-Hn2[-19190.7 - 2426.57 6^/4 - 317.732 ^''/16 - 1547.38 5^/64] 
-|-n3[62169.5 + 2521.29 6^/4 + 3470.28 6'^/16] 

-|-n'^[-91005.1 + 3984.82 S^/4 - 9148.6 6^/16]] (Al) 

and the vector field in MeV by 

I]°(r, n, d) = 11° (T, n, -5) ^ n[3462.24 + 946.705 5/2 - 0.334508 5^/4] 

+n2hll312.4- 6246.21 6/2 - 6353.53 6^/4 ~ 0.099478 6^/8] 

[20806.1 -I- 18717.6 6/2 + 29298. 6^/4 - 0.490543 6^/8] 
+n*[352.371 - 24887.2 6/2 - 39807.4 (5V4 - 0.346218(5^8]. (A2) 

The vector field is nearly independent of temperature, the scalar field has a weak temperature dependence. 



